Statistical geometry in scalar turbulence 
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A general link between geometry and intermittency in passive scalar turbulence is established. 
Intermittency is qualitatively traced back to events where tracer particles stay for anomalousy long 
times in degenerate geometries characterized by strong clustering. The quantitative counterpart is 
the existence of special functions of particle configurations which are statistically invariant under 
the flow. These are the statistical integrals of motion controlling the scalar statistics at small scales 
and responsible for the breaking of scale invariance associated to intermittency. 



PACS number(s) : 47.10. +g, 47.27.-i, 05.40. +j 

Scalar fields transported by turbulent flow occur in 
many physical situations, ranging from the dynamics of 
the atmosphere and the ocean to chemical engineering 
(see, e.g., Ref. Q). Specific examples are provided by 
pollutant density, temperature or humidity fields and 
the concentration of chemical or biological species. The 
advection-diffusion equation governing the transport of 
the scalar field 8 is : 



d t 9(r, t) + (v V) 6(r, t) = K A6(r, t), 



(1) 



where v(r,t) is the incompressible advecting flow and 
k is the molecular diffusivity. Two broad cases are dis- 
tinguished : active scalars, where v depends on 9, e.g. 
by an explicit relation v = v (6), and passive scalars, 
where the statistics of v is independent of 9. Here, we 
shall be concerned with the latter, although we conjec- 
ture that the physical mechanisms presented in the fol- 
lowing are quite general and relevant also for the active 
cases. The Fokker-Planck equation (|l|) is associated to 
the Lagrangian dynamics of tracer particles whose posi- 
tion p(t) obeys dp(t) = v(p(t),t)dt + \/2k df3(t), where 
(3(t) is the isotropic Brownian motion The equation 
(|l|) governs the evolution of the probability density of 
particles at position r and time t. 

Scalar turbulence is typically generated by maintaining 
a mean scalar gradient (9) = gr, e.g. by heating/cooling 
devices in temperature field experiments. The notation 
(•) denotes the average with respect to the velocity statis- 
tics, which is in principle arbitrary. We shall be inter- 
ested in flows with correlations having a nontrivial power 
law behavior in the inertial range of scales r <C L, where 
L is the velocity correlation length. Examples are pro- 
vided by the two and three dimensional Navier-Stokes 
turbulent flow (see, e.g., Ref. ||). A very robust fea- 
ture of scalar turbulence is its strong intermittency : rare 
strong events (such as the sharp cliffs observed in Fig. 1) 
dominate the scalar statistical properties. More quanti- 
tatively, intermittency reflects in the anomalous scaling 
of the correlations. In the inertial range, the scalar struc- 
ture functions S n (r) = ((9(r, t) — 9(0, t)) n ) take the form 



with a nonvanishing value of C» lm — Cn- Here, Cn™ is the 
value predicted by simple mean field dimensional argu- 
ments, e.g. of the Kolmogorov 1941 type. The positive 
value of the anomalous correction C* m ~ Cn reflects the 
breaking of scale-invariance : the scale L explicitly ap- 
pears in the inertial range expressions of scalar correla- 
tions, even though r <C L. The phenomenon of intermit- 
tency is quite generic for scalar turbulence and indepen- 
dent of the specific choice of v, including for Gaussian 
velocity fields (see, e.g., Ref. 0). Past research on inter- 
mittency has mostly concentrated on phenomenological 
models making vague contacts with the dynamics and 
based uniquely on scaling exponents and structure func- 
tions. Angular dependencies and geometrical informa- 
tions about the correlations were discarded. It is shown 
here that the whole structure of multi-point correlations 
is in fact needed to achieve a real dynamical understand- 
ing based on the geometry of the figures identified by 
tracer particles. Deviations of exponents from naives di- 
mensional values turn out to be just by-products of the 
figure geometry nontrivial evolution under Lagrangian 
dynamics. While the figure size grows according to sim- 
ple dimensional arguments, the evolution of their shape 
is a more delicate issue. There are in particular some spe- 
cial functions of the particle positions whose average with 
respect to the Lagrangian dynamics remains constant in 
time. That is due to a delicate compensation between the 
growth due to the figure size and the depletion associated 
to the figure shape. Those statistically conserved func- 
tions are dominating the behavior in the inertial range 
and controlling the scalar field intermittency 

Specifically, the velocity field v considered here is a 
two-dimensional turbulent flow generated by an inverse 



energy cascade process 
both for experiments [ 



S„(r)«r<» m f- 



(2) 



This is a situation of interest 
and in the atmosphere 
The flow realizes the type of turbulence theoretically pos- 
tulated by A.N. Kolmogorov in 1941 : it is isotropic, it 
has a constant energy flux (but upscale) and it is scale- 
invariant with scaling exponent 1/3 A property of 
interest to us is that the velocity is not intermittent. All 
nontrivial scaling properties of the scalar field presented 
in the following are therefore entirely due to the 
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FIG. 1. A snapshot of the scalar field obtained by numer- 
ical integration of (1) with advection by a two-dimensional 
turbulent flow generated by an inverse energy cascade. The 
scalar turbulence is maintained by a fixed temperature gradi- 
ent (0) — g ■ r, with g oriented from right to left. 

advection-diffusion equation (|l|) and not mere footprints 
of the velocity field. Details on the integration proce- 
dure used for the numerical simulations of (Q) with a 
mean gradient g can be found in Ref. jlO|. The single- 
time scalar statistics at the stationary state is defined by 
the n-point correlations C n (r,t) — (0(ri, t) ■ ■ ■ 6(r n , t)), 
where r denotes the set n, . . . , r„. For spatially homoge- 
neous situations, C n is invariant under translations and 
it depends only on 2n — 2 degrees of freedom, associated 
to the separation vectors among the n points. In the 
inertial range, the velocity scale invariance is expected 
to reflect in scalar correlations C n (r) behaving as power 
laws with respect to global size variables, such as e.g. the 
gyration radius of the set r. 

Let us consider for simplicity the third order case n — 
3. All the following arguments are easily generalized to 
higher order correlations. The correlation function C3 
depends on the size, the orientation and the shape of the 
triangle defined by the three points n, r 2 and r%. The 
global size variable can be defined as R 2 — (r\ 2 + r\\ + 
r| 1 )/3, where r,y is the distance between the i-ih and the 
j-th particle. As shown in Fig. 2, in the inertial range of 
scales C3 depends on R as a power law with the exponent 
(3 = 1.25. The hallmark of intermittency is in the fact 
that £3 is smaller than the dimensional prediction 5/3 
(see Ref. As for the shape and the orientation of the 
triangle, we shall use the same Euler parametrization as 
in Refs. PO]. Defining p 1 = (n - r 2 ) /V2 and 



FIG. 2. The dependence of the third order correlation func- 
tion C3 with respect to the size of the triangle R. The straight 
line is the power law behavior R 125 . 

p 2 = (ri + r 2 — 2r 3 ) /V&, the shape of the triangle is 
controlled by the two variables 



X 



1/2 tan" 
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w = 2 



I Pi x P2I 
R 



(3) 



Some of the shapes associated to different values of \ 
and w are shown in Fig. 2. The global orientation of the 
triangle with respect to the mean gradient direction g is 
defined by the angle ip. It is convenient to decompose 
C3 on the orthogonal basis made of cos (£ip) and sin (£<p) . 
Reversing the coordinates with respect to an axis parallel 
or orthogonal to g statistically leaves the 9 field invariant 
or inverts its sign, respectively. In the projection of C3, 
the angular momentum £ should therefore be odd and 
sine functions are absent. Furthermore, the dominant 
contribution at the small scales is the one having the 
lowest angular momentum (see Ref. |l3| for the case of 
Navier-Stokes turbulence). The correlation function C3 
takes then the form 



C 3 (r) = R Ca fix, to) cos ip + 



(4) 



where the dots stand for subdominant higher-order har- 
monics of the form cos(2£ + l)tp. The invariance under 
arbitrary permutations of the three vertices of the trian- 
gle allows to reduce the phase space to — tt/6 < x < tt/6, 
< w < 1 and the function / in (^) is antiperiodic in \ 
with period 7r/3 P,Q. The measured dependence of / 
on the shape coordinates \ and w is shown in Fig. 3. The 
maximum of / is realized at \ — 0, w = 0, where two 
of the three particles are stuck together. For equilateral 
triangles (w = 1) or for the "dumbbell" configuration 
X = tt/6, w = 0, the symmetries enforce / = 0. 

Let us establish the connection with the Lagrangian 
dynamics. That is done by replacing the Eulerian vari- 
ables r in the argument of C3 by their Lagrangian evo- 
lutions p(t), with p(0) = r. The resulting object is a 
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FIG. 3. Contour lines in the \ — w plane of the third order 
shape function / appearing in (5). 

stochastic function whose average with respect to the La- 
grangian trajectory statistics is denoted by (•)/:■ More 
generally, for a generic function <fi(r) of the n points Ti 
we can define its Lagrangian average as 



{<i>{t))c = J 4>{p)Pn (t,p\0,r) dp. 



(5) 



Here, the n-particle propagator P n (t , p|0 , r) denotes the 
probability that, being in r at time 0, the n particles are 
in p at time t. In a turbulent Kolmogorov flow, distances 

typically grow with time as |i| 3 ' 2 , whose special instance 
is the celebrated Richardson law (r 2 2 ),c(i) °c |i| 3 for the 
distance r\i between two particles. For generic functions 
0, homogeneous of positive degree a that is 4>(Xr) = 
A cr 0(r), the Lagrangian average is therefore expected to 
grow as |£| 3<T / 2 . 

Intermittency is dynamically originated by blatant vi- 
olations of expectations a la Richardson. The anomalous 
part of correlation functions has in fact a constant La- 
grangian average, as clearly demonstrated in Fig. 4 for 
the third order case. Those statistically preserved func- 
tions are the statistical integrals of motion responsible 
for the breaking of scale invariance associated to inter- 
mittency. The other important point is that their con- 
stancy is tightly related to the geometry of the figures 
identified by the Lagrangian particles. Fig. 4 indicates 
indeed that the size factor R^ 3 grows as |i| 3l » 3 / 2 . The 
Lagrangian average of C3 remaining constant, the shape 
part f{x> w ) cos V m © must compensate for the growth 
of the figure size. As shown in Fig. 3, the function / de- 
creases going from degenerate triangles with two of the 
vertices close to each other to triangles with aspect ratios 
of order unity. The geometrical meaning of anomalous 
scaling laws is then quite clear: the smaller is £3, the 
slower is the compensation needed from the shape factor 
and the longer degenerate triangle configurations persist. 
In other words, "stronger intermittency = particles stay- 
ing longer close to each other" . 



FIG. 4. The Lagrangian average of the correlation function 
{C r i)c- For comparison, it is also shown the evolution of the 
average {R c,3 )l ~ |£| 5< " 3 , that obeys the dimensional scaling 
law. 

Systematic support to the previous physical ideas can 
be provided in the special case where the advecting ve- 
locity in ([!]) has a short correlation time, the so-called 
Kraichnan model |l4j]. The assumption is of course 
far from realistic, but it leads to the peculiar property 
that the n-particle propagators P n in (0) obey closed 
Fokker-Planck equations |l5| (see also Ref. Q). The sta- 
tistically preserved functions are now identified as zero 
modes of the n-particle Fokker-Planck operator [|l6] [is) . 
Their anomalous scaling behavior could be calculated in 
some perturbative limits ||l6|-|l8|| or measured numerically 
|l9| , p0| . Zero modes enter the P n 's via the asymptotic ex- 
pansion H] (see also Ref. p3) : 



Pjt,p 



0,Ar)=^ ^M*,£), 



(6) 



valid for small A's. Zero modes are the q = terms 
in (^) and they are ordered according to their scaling 
dimension by the index i. Higher q's identify the so- 
called slow modes, whose Lagrangian average is growing 
as an integer power of time, although with an exponent 
smaller than the dimensional one 3<7j jg /2. 

A simple example of slow mode for our inverse cascade 
velocity field is provided for n = 2. The Lagrangian 
average of (g ■ 7*12) is preserved as its time derivative is 
proportional to ((i>i — t^)) = 0. The first slow mode 
associated to it is given by the function (g ■ r^) r\(^ . 
Fig. 5 shows indeed that its Lagrangian average grows as 
|t|, much slower at large times than the dimensional law 

|i| 5 / 2 . 

What is the degree of generality of the Lagrangian 
preservation mechanism for intermittency and the expan- 
sion (Q) ? The crucial property of the Kraichnan velocity 
field is its short correlation time (this ensures that the 
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FIG. 5. The Lagrangian average of the anisotropic slow 
mode (g ■ ru) r 

12 time. 

Lagrangian trajectories are Markovian). For the inverse 
energy cascade flow considered here this property is lost 
as the correlation time is finite. The numerical results 
presented here give therefore a strong indication that the 
basic mechanisms for scalar intermittency are quite ro- 
bust and generic. Changing the statistics of the flow 
affects only quantitative details, such as the numerical 
value of the exponents or the precise shape of the La- 
grangian preserved functions. The expansion (|^) is analo- 
gous to those encountered in many-body statistical prob- 
lems p|| and controls the small-scale behavior of scalar 
correlations. Indeed, it follows from ([!]) that the scalar 
is preserved along the Lagrangian trajectories. Scalar 
correlations can be then expressed as 

C„(r,t) = J P n (~t,p\0,r) (g- Pl )...(g- Pn )dp. (7) 

Inserting @ into this expression it is evident that the be- 
havior at small scales is a superposition of the functions 
4>i,q- Two general remarks on scalar turbulence follow. 
First, changing the initial condition g ■ r (or the injection 
mechanism in the forced case) modifies the constants but 
not the scaling exponents in the correlation functions. 
This naturally explains the universality properties of 
scalar turbulence observed in Ref. |l(J. Second, all fewer- 
particle modes (f>i, q , i.e. those depending on m < n vari- 
ables, appear in the expression for C n . Their Lagrangian 
averages with respect to the evolution of the m particles 
or the whole set of n particles are indeed trivially coin- 
ciding. Note however that structure functions satisfy the 
trivial identity S n (r) = f ■■■ f d n ■ ■ ■ d rn C n (r) dr. All 
m-particle modes will therefore drop out from the n-th 
order structure function. This clearly illustrates a point 
previously mentioned : the apparent simplification of the 
single distance r left in structure functions is purely illu- 
sory as their anomalous scaling laws are still dynamically 
associated to n-particle geometries. 



In conclusion, we have identified the integrals of motion 
responsible for the failure of dimensional analysis and in- 
termittency in scalar turbulence. The anomalous parts of 
correlation functions are preserved under the Lagrangian 
dynamics. That preservation is a geometrical effect asso- 
ciated to figures identified by scalar particles persisting 
in degenerate shapes. Particles tend indeed to remain 
close to each other for unexpectedly long times, pointing 
to applications for the reaction rate enhancement in the 
transport of chemically reactive species. 
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